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i— 1 ; ABSTRACT 

5-H , In a coalescing neutron star-neutron star or neutron star-black hole binary, oscillation 

modes of the neutron star can be resonantly excited by the companion during the 
final minutes of the binary inspiral, when the orbital frequency sweeps up from a few 
Hertz to a few thousand Hertz. The resulting resonant energy transfer between the 
orbit and the neutron star speeds up or slows down the inspiral, depending on whether 
^v^j \ the resonant mode has positive or negative energy, and induces a phase change in the 

■ emitted gravitational waves from the binary. While only g-modes can be excited for a 

' nonrotating neutron star, f-modes and r-modes can also be excited when the neutron 

t-H , star is spinning. A tidal resonance, designated by the index (J k, m) ({jk} specifies the 

t-H ■ angular order of the mode, as in the spherical harmonic Yjk), occurs when the mode 

| frequency equals m times the orbital frequency. For the f-mode resonance to occur 

before coalescence, the neutron star must have rapid rotation, with spin frequency 
v s <; 710 Hz for the (22, 2)-resonance and v s <; 570 Hz for the (33, 3)-resonance 
(assuming canonical neutron star mass, M = 1.4 Mq, and radius, R = 10 km; however, 
<** j ■ for R = 15 km, these critical spin frequencies become 330 Hz and 260 Hz, respectively). 

Qh Although current observations suggest that such high rotation rates may be unlikely 

q for coalescing binary neutron stars, these rates are physically allowed. Because of their 

5^ ■ strong tidal coupling, the f-mode resonances induce a large change in the number 

| of orbital cycles of coalescence, AA or b, with the maximum AA or b ~ 10 — 1000 for 

the (22, 2)-resonance and AA or b ~ 1 for the (33, 3)-resonance. Such f-mode resonant 
effects, if present, must be included in constructing the templates of waveforms used in 
. £^ | searching for gravitational wave signals. Higher order f-mode resonances can occur at 

. slower rotation rates, but the induced orbital change is much smaller (AiV or b 0.1). 

' For the dominant g-mode (22, 2)-resonance, even modest rotation (v s <, 100 Hz) can 

enhance the resonant effect on the orbit by shifting the resonance to a smaller orbital 
frequency. However, because of the weak coupling between the g-mode and the tidal 
potential, A./V or b lies in the range 10~ 3 — 10~ 2 (depending strongly on the neutron star 
equation of state) and is probably negligible for the purpose of detecting gravitational 
waves. Resonant excitations of r-modes require misaligned spin-orbit inclinations, and 
the dominant resonances correspond to octopolar excitations of the j = k = 2 mode, 
with (jk,m) — (22,3) and (22,1). Since the tidal coupling of the r-mode depends 
strongly on rotation rate, AJV orb ^ 10 _2 (i?/10 km) 10 (Af/1.4 Mq) -20 / 3 is negligible for 
canonical neutron star parameters but can be appreciable if the neutron star radius 
is larger. 

Key words: binaries: close - gravitation - hydrodynamics - stars: neutron - stars: 
oscillation - stars: rotation 
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1 INTRODUCTION 

The new generation of gravitational wave detectors such 
as the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO; Abramovici et al. 1992) and its French-Italian 



counterpart VIRGO (Bradaschia et al. 1990) are designed 
to observe gravitational radiation in the frequency range 
from ~ 10-1000 Hz. Coalescing neutron star-neutron star 
(NS-NS) and neutron star-black hole (NS-BH) binaries are 
promising sources for such radiation, with the last few min- 
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utes of the binary inspiral producing gravitational waves 
which sweep upward through this frequency range (e.g., 
Thorne 1987, Cutler et al. 1993, Thorne 1998). Due to the 
low signal-to-noise ratios, a good understanding of the mech- 
anisms that generate and affect the gravitational waveform 
is required to create accurate theoretical templates that can 
be used to extract the signals from the noise. 

The binary inspiral and the resulting waveform can be 
described, to leading order, by Newtonian dynamics of two 
point masses, together with the lowest order dissipative ef- 
fect corresponding to the emission of gravitational radiation 
via the quadrupole formula. To more accurately construct 
the theoretical templates, general relativistic effects must 
be taken into account by using post-Newtonian expansions 
(e.g., Blanchet et al. 1995) or techniques of numerical rela- 
tivity (e.g., Teukolsky 1998; Brady et al. 1998). Other cor- 
rections come from hydrodynamical effects due to the finite 
size of neutron stars (Kochanek 1992; Bildsten & Cutler 
1992; Lai et al. 1994). In particular, tidal interactions can 
have significant effects on the final stage of binary inspiral 
by destabilizing the orbit and accelerating the coalescence at 
small orbital radii (e.g., Lai et al. 1994; Lai & Shapiro 1995; 
Rasio & Shapiro 1995; Lai & Wiseman 1996; Baumgarte 
et al. 1997; New & Tohline 1997; Uryu & Eriguchi 1998; 
see also Shibata et al. 1993; Davies et al. 1994; Zhuge et 
al. 1996; Ruffert et al. 1997 for numerical simulations on 
binary coalescence). 

The early stage of the inspiral, with the gravitational 
wave frequencies in the range of 10 Hz to a few hundred Hz, 
is particularly important for detecting the wave signal and 
for extracting binary parameters from the waveform using 
the matched filtering technique (e.g., Cutler et al. 1993; Cut- 
ler & Flanagan 1994) . It is usually thought that tidal effects 
are completely negligible for this "low- frequency" regime. 
This is indeed the case for the "quasi-equilibrium" tides as 
the tidal interaction potential scales as D~ 6 (where D is the 
orbital separation; see Lai et al. 1994 for analytic expres- 
sions of the orbital phase error induced by static tides). The 
situation is more complicated for the resonant tides: As two 
compact objects inspiral, the orbit can momentarily come 
in resonance with the normal oscillation modes of the NS. 
By drawing energy from the orbital motion and resonantly 
exciting the modes, the rate of inspiral is modified, giving 
rise to a phase error in the gravitational waveform. This 
situation has been studied (Reisenegger & Goldreich 1994; 
Lai 1994) in the case of nonrotating neutron stars where 
the only modes that can be resonantly excited are g-modes 
(with typical mode frequencies ^ 100 Hz). It was found that 
the effect is negligible because the coupling between the g- 
mode and the tidal potential is weak. However, the NS can 
be rotating, and consequently, these results need to be re- 
examined. As we show in this paper, even modest rotation 
of the NS can enhance the resonant g-mode excitation by 
shifting the resonance to a lower frequency (where the or- 
bital decay is slower and resonance is stronger). In addition, 
f-modes, which possess strong tidal coupling, can have their 
normally high frequencies in the nonrotating case "dragged" 
or reduced by retrograde rotation to sufficiently low frequen- 
cies in the inertial frame, such that resonant excitations can 
occur well before binary merger. The NS spin also gives rise 
to a new class of modes which are absent in the nonrotating 



case: R-modes, which can have low frequencies and may also 
undergo resonant excitations. 

The purpose of this paper is to study the resonant ex- 
citations of various modes of spinning NSs in inspiraling 
binary systems and to examine their effects on the orbital 
decay rate (particularly the phase error in the gravitational 
waveform). Section 2 develops the general formalism of res- 
onant tides in coalescing NS binaries. In section 3 we apply 
the analysis to f-modes, g-modes, and r-modes when the 
spin axis is aligned or anti-aligned with the orbital angu- 
lar momentum axis. Section 4 describes the modification to 
the results of the aligned (anti-aligned) case when the spin 
and orbital axes are misaligned by an inclination angle (3. In 
section 5, we summarize our results and discuss their impli- 
cations for gravitational wave detection. 



2 RESONANT MODE EXCITATION 
2.1 Basic equations 

Consider a NS with mass Al, radius R, and spin $~2 S in orbit 
with a companion Af', which can be another NS or a BH. 
We treat M' as a point mass. The orbital frequency vector 
is fiorb- We define a coordinate system (x,y,z) centered on 
M with the z-axis along f2 or b, such that f2 or b = O rbZ (with 
f^orb > 0). Define another coordinate system (x',y',z') with 
the z'-axis along tt 3 , such that tt 3 = £l a z' (with Q s > 0) and 
the j/'-axis in the orbital plane; these are the body axes. Let 
the angle between the z-axis and z'-axis be /3 (the spin-orbit 
inclination angle) and the angle between the y-axis and y'- 
axis be a. Thus the (x' , y' , z') frame is related to the (x, y, z) 
frame by Euler angle (a,/3, 7 = 0). 

The gravitational potential produced by M' can be ex- 
panded in terms of spherical harmonics: 



U(r,t) = —. 



GM' 



-D(t)| 

-gm'y: 



Wlr 



Mt)] 



(2.1) 



where r = (r, 8, <f>) is the position vector of a fluid element in 
star M (r, 6, 4> are the spherical coordinates with respect to 
the z-axis), D(t) = [D(t), n/2, 3?(t)] is the position vector of 
the point mass, M', and Wi m is defined in Press & Teukolsky 
(1977) as 



= (-) 



(i + m)/2 



47T 

27+1 



(l + m)l(l-m)\ 



1/2 



I + m 



I — m 



(2.2) 



Here the symbol (— ) p is zero if p is not an integer. For 
example, W 2 ±2 = (37r/10) 1/2 and W 3 ±a = (5tt/28) 1/2 . In 
equation (2.1), the I — and 1 = 1 terms can be dropped 
since they are not relevant for tidal deformation. The tidal 
potential is specified by the index {Irn} (relative to the or- 
bital angular momentum axis z) and is nonzero only when 
(l+m) is even, reflecting the symmetry of the tidal potential. 

It is more convenient to define oscillation modes rel- 
ative to the spin axis z', and we will need to express the 
tidal potential in terms of Yi m i (9', <f>'), the spherical har- 
monic function relative to z'. The function Yi m {6, <f>) is re- 
lated to Y lm ,(9',<j>') by 
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Y lm (6,$) = ^^» m ( a ,i3)7 im ,(ff'^') 



(2.3) 



where the Wigner ©-function is given by (e.g., Wybourne 
1974) 



V^ m (a, P) = e m ' a [(I + m)\ (I - m)\ (I + m')\ (I - m')\ 



1/2 



E 



«+™'-p (cosf) 2p " m " m '(smf) 2; - 2p+m+m ' 
p\ (I + m — p)\ (I + m' — p)\ (p — m — m')\ 



(-1) 



(2.4) 



Thus {lm 1 } is the tidal potential index relative to the spin 
axis. 

The linear perturbation of the tidal potential on M is 
specified by the Lagrangian displacement, £(r,t), of a fluid 
element from its unperturbed position. In the inertial frame, 
the equation of motion takes the form 



d 2 £ d£ 



(2.5) 



where v = tt 3 x r is the unperturbed fluid velocity and C is 
a self-adjoint operator (a function of the pressure and grav- 
ity perturbations) acting on £ (Friedman & Schutz 1978a). 
Our analysis is accomplished through an eigenmode decom- 
position of the Lagrangian displacement: £ is taken to be 
the linear superposition of eigenmodes, £ a (r,t), namely 



(2.6) 



We shall specify the normal mode by the index a = {njk} 
with respect to the body axes: n is the number of radial 
nodes, j gives the number of nodes in the S'-eigenfunction, 
and k is the azimuthal index (about the z'-axis). In the 
following we shall suppress the radial index n so that a — 
{jk}. In the absence of external perturbation, the eigenmode 
behaves as 



|* Q (r,t)=|* Q (r) e ^ t oc e l ^ t+l ^', 



(2.7) 



where ui a is the a-mode frequency in the inertial frame. The 
mode eigenfunction satisfies the equation 

- uUc + 2iu a (v ■ V% + C • £» = 0, (2.8) 

and we shall adopt the normalization 



d X p ^ a • &aa' • 



(2.9) 



Substituting (2.6) into (2.5), using (2.7) and (2.9), and 



noting that (2.8) gives a relation for C ■ £, the dynamical 



equation for the mode amplitude is 



d a — 2iT a a a + (u) a — 2T a ui a )a a 

G 
D 



E GM' (l) _ im q> 



(2.10) 



Irara' 



where 



i d x pC a - (v ■ V)£ Q 



kQ 3 +i I d?x pCa- (n s X £ a ) 



(2.11) 



(this is called B a in Lai 1997) and Q a ^ m ' is the tidal cou- 



pling coefficient or overlap integral defined by 

Qa.lm' = j d\' p&-V[r l Y lm , 

r d 3 x'Sp* a [r l Y lm ,(8',^')]. (2.12) 

Here Sp a — —V • (p£, a ) is the Eulerian density perturbation. 
Clearly, Q a ,im' is nonzero only when m' = k, thus 



Qa,lm' Qjk,lk^m'k~ 



2.2 Resonant Mode Energy 



(2.13) 



The right-hand-side of equation (2.10) contains a sum 
of driving terms with time dependence proportional to 
e"'™*' 1 '. For a given mode a = {jk}, resonance occurs when 



\ui a \ = |m|57 or b, (Resonance Condition) 



(2.14) 



provided there exists at least one I (obviously I > \m\, I > 
\k\) which satisfies 

I + m — even (Symmetry of the Tidal Potential) (2.15) 



and 

r D k l) m Qjk,ik + 0. (Nonzero Tidal Coupling) 



(2.16) 



For Qjk,ik to be nonzero, the projected tidal potential [oc 
Yik{9' , (j>')\ must have the same ^'-symmetry, with respect 
to the equator, as the mode. This gives the requirement 
that I + j = even (for spherical modes) or I + j = odd (for 
r- modes) (see §3). 

We consider now a specific resonance. Hereafter we shall 
adopt the convention^] k > and m > 0. Thus a resonance 
is designated by the indices (a,m) = (jk,m), which indi- 
cates that the mode a — {jk} is excited by tidal potential 
of azimuthal order m (with angular dependence e ±!m *). In 
general, there can be many Vs that satisfies the conditions 
( [2.14| )-( |2.1(f )p], but clearly the smallest I which satisfies these 
conditions will provide the dominant contribution. When ap- 
propriate, we shall also label a resonance using the notation 
(jk, lm), with / specifying the angular order of the dominant 
tidal field that contributes the resonance (jk,m). 

In the special case when the spin and orbit are aligned 
or anti-aligned, one must have m — k, i.e., each mode can 
experience resonance only once. But in the general case of 
arbitrary spin-orbit inclination angle, f3, for a given mode 
{jk}, there may b e infi nite number of m which satisfies the 
conditions (2.14)-(2.16); in other words, the mode can be 
resonantly excited many times at different orbital frequen- 
cies. Thus, strictly speaking, the energy transfer from the 
orbit to the mode at a given resonance (jk, m) depends on 
the amplitude of the mode induced by previous resonances 
(jk, m + 2), (jk, m + 4), • • • 0. However, as we shall see be- 
low, in almost all cases of interest, the energy transfer at the 



k and m can be positive 



* Clearly, in equations (2.7) and ( [2.1 
or negative. 

T The exception is when the star is nonrotating, in which case 
m = k and I = j. 

t Consider a spheroidal (f or g) mode, a = {jk}. For the (jk,m)- 
reso nance to be in effect, it is necessary for I + m = even [see 
eq. (3.15)] and l + j = even (otherwise Qjk,ik = 0)- Now in order 
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(jk, m)-resonance, AEjk, m , is much greater than the energy 
transfer at the earlier resonances, i.e., AEjk, m 2> AEjk,m+2- 
Therefore, throughout this paper, we shall treat different 
resonances of the same mode as being independent of each 
other. This corresponds to assuming that before each reso- 
nance, the mode amplitude is essentially zero. 

To calculate the energy transfer to a given mode at a 
particular m-resonance, A Ejk.m , we only need to keep the 
resonant driving terms in (2.10), namely 



2T a uj a )a a 



(0 " a , h e ~ im * 



e im *. (2.17) 



l 



Sinc e a a oc e 1 ^" in the absence of external forcing [see 
eq. (2.7)], it is natural to define the mode amplitude c Q via 



(2.18) 



and we expect c a to be independent of time after the reso- 



nance Q Equation (2.17) then becomes 
c Q + 2i(w Q - T a )c a 

(\~* GM' (n \ _j m $_j 



(2.19) 



Only one of the terms on the RHS of (2.1£) can induce the 
resonance, which occurs when the argument of the exponen- 
tial of either the first term or the second term goes through 
zero. The term that induces resonance depends on the sign 
of bio.. 



2.2.1 General Properties of the Mode Frequencies 

Before proceeding with our calculation, let us briefly discuss 
various possibilities for the mode frequency. In general, for a 
given mode index a = {jk} (recall our sign convention k > 
0, and note that we have suppressed the radial mode index 
n) , there are two possible solutions of opposite signs for wi^ , 
the mode frequency in the frame corotating with the star, 
which correspond to two different dire ctions of propagation 
for the wave pattern. Recall from (2.7) that our convention 



for free oscillations is e 



iuj a t-\-ik4> 



relative to the spin axis, 



where fl a — Q s z' and Q, 3 > 0. The mode frequency in the 
inertial frame, uj a , is related to oj£'' by 



(2.20) 



1 



for the (jk, m + 1) to be in effect, it is necessary for I + m - 
even, which is incompatible with I + j = even. 
§ Before and after a resonance, there is a small, non-resonant 
contribution to a a which is proportional to e ± Im *(*). This con- 
tribution can be neglected compared to the resonance-induced 
amplitude. 



(r) 

The first solution, with ui a > 0, corresponds to the 
mode which propagates retrograde with respect to the spin 
(in the rotating frame). Rotation drags the mode frequency 
in the inertial frame to small values; thus the resonance is 
shifted to smaller orbital frequency. With increasing Q a , the 
mode frequency ui a may become negative, corresponding to 
wave propagation in the same direction as the spin (in the 
inertial frame) . This change in sign of uj a signals the onset 
of the so-called Chandrasekhar-Friedman-Schutz (CFS) in- 
stability (Chandrasekhar 1970; Friedman & Schutz 1978b) 
and can occur for both f- modes (e.g., Lindblom 1995) and 
g- modes (Lai 1998). At this point the mode energy in the 
inertial frame changes from positive (stable mode) to neg- 
ative (unstable mode)^. Note that when the resonant exci- 
tation of a positive-energy mode occurs, the orbital energy 
is transferred to the mode and the orbit decays faster; when 
a negative-energy mode is excited, the orbit actually gains 
energy and the orbital decay is slowed. 

The second solution, with LOa ' 1 < 0, propagates pro- 
grade with respect to the spin (in the rotating frame). Ro- 
tation causes u) a to become even more negative. This mode 
always has positive energy and is stable. The resonance is 
shifted to even higher orbital frequency. Thus for f-modes, 
the uJa ^ < modes do not experience any resonance. For 
g- modes, |w a | can still lie in the appropriate range for reso- 
nance. 

(r) 

The r-modes are a special case for which u a > but 
uj a < 0, meaning that in the rotating frame, the r-mode wave 
pattern always propagates opposite to the spin, while in the 
inertial frame it propagates in the same direction as the spin. 
This implies that the r-modes are always CFS unstable in 
the absence of fluid viscosity (see §3.3 below). 



2.2.2 Mode with uj a > (Spin-Retrograde Mode) 

For a mode with oj a > 0, which corresponds to wave 
propagation opposite to the rotation, resonance occurs at 
m"l> = mSl or b = uia, where Q Q ib is given by 



2 _G(M + M') 

* 'orb 



D'- 



(2.21) 



neglecting tidal and post-Newtonian cor rectio ns. In this 
case, only the second term on the RHS of ( 2.19 ) can induce 
resonance. Thus we have 



c a + 2i(u a — T a 



E 



^WW2?« m Q aiI J e*™*-***. (2.22) 



Integrating over time around the resonance and noting that 
c Q is almost constant after the resonance (and c Q = before 
the resonance since the mode amplitude induced by previous 
resonances is negligible; see the discussion at the beginning 
of §2), we find 

1 

c a (after resonance) 



2i(u; a — T a ) 



^ Since the growth time of the instability can be very long, and 
since viscosity can suppress the CFS instability for some modes 
(e.g., Lindblom 1995), we shall still consider such unstable modes 
in this paper. 
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(2.23) 



Resonance occurs when the argument of the exponential is 
equal to zero. Doing a Taylor expansion of the argument 
of the exponential around the resonance epoch and using 
m$ = u) a , we have 

m$ — uj a t — m&t + mi$t 2 — cu a t — ^Cl OI bt 2 . (2.24) 
D varies slowly with time, so that 



dt ■ 



1 



2?r 



1/2 



(2.25) 



Therefore the mode amplitude after the («, m)-resonance is 



a a {t) 



2i(uj a ~T a ) 



2tt 



1/2 



, (2.26) 



where D a is the orbital radius where resonance occurs and 
f2 or b is evaluated at D a . The energy of the mode in the 
inertial frame is given by Friedman & Schutz (1978a) as 



+ r • c ■ i 



(2.27) 



Using eqs. (2.6)-( p.9[ ), this gives the mode energy associated 
with the (a, m)-resonance as 

= _> 

= 2ui a (ui a - T a )\a a \ 2 



AE at „ 



\d a \ 2 + -(wo — 2T a Lo a )\a a \ 2 



0J a — T a 



mtt olh J a \ 



(2.28) 



The factor of two in the first line of ( 2.2S ) arises because, 
for each k > and uj a > resonant mode, there is also 
the identical k < and uj a < resonant mode (recall our 
sign convention m > and k > 0). Hereafter we shall defin e 
Qa,ik in units such that G = M = R = 1. Equation ( 2.2S ) 
then becomes 



AE a , 



(GM_\ ( GM' 2 



1 Q 



;+i 



Here we have defined 

rji (r) 



d\pC a ■ (n, X f a ) 



(2.29) 



(2.30) 



Since the effect of the tidal interaction on the orbit is small, 
we use the quadrupole formula for the gravitational radi- 
ation of two point masses to obtain the orbital decay rate 
(e.g., Shapiro & Teukolsky 1983) 



6 -_ 3 ^o n_ 64G 3 Af 3 g(l + g) 

"orb- 2I) "orb, u- 5c5 D3 
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(2.31) 



where q = M'/M is the mass ratio. Equation (2.29) can be 



expressed in terms of index (a, m) = (jk, m) and quantities 
associated with the resonant mode, giving 



96 V R I \GM 



1 /q> a V 

(1 +fl)(2I+l)/6 V m J 



2(!-l)/3 



where we have defined dimensionless quantities 

53 \ 1/2 



Y gm y 



(- 

\GM 



3 \ 1/2 



(2.32) 

(2.33) 
(2.34) 
(2.35) 



In this case, the mode energy is always positive (see §2.2.1). 



2.2.3 Mode with uj a < (Spin-Prograde Mode) 

For a mode with u a < 0, corresponding to wave propaga- 
tion in the same direction as the spin, the resonance occurs 
at mfl or b = — u) a , and only the first term on the RHS of 
cq. ( ]2 . 1 E| ) can induce resonance. Carrying out a similar cal- 
culation as in the u) a > case, we find 



AE a .m — 



96 V R J \GM 



5tt (GM^\ (Rc 2 Y' /2 g 



(l + 9 )(«+i)/a 



2(i-l)/3 



W lm Dl 



m 



(2.36) 



Note that the mode energy can be either positive or negat ive, 
depending on the sign of ui a e a [see §2.2.1 and eq. (2.28)]. 



Apart f rom the sign in uj a , the only difference between 



and ( 2.36) is the factor T> k 



(0 



and V { l ] m . This differ- 



ence reflects the "preferred" spin-orbit orientation for reso- 
nant excitation: In order to have resonance, the wave pattern 
of the mode in the inertial frame must propagate in the same 
direction as the orbital motion of the external mass, M' . For 
uj a > 0, the wave propagates opposite to the spin, thus the 
resonance will be strongest for retrograde spin-orbit incli- 
nation (/3 > 90°); For uj a < 0, the wave propagates in the 
same sense as the spin, thus we expect the resonance to be 
strongest for prograde spin-orbit inclination (/3 < 90°). 

2.3 Change in Orbital Phase 

We now calculate the orbital phase error due to the resonant 
energy transfer between the binary orbit and the NS oscilla- 
tion mode. The number of orbital cycles, N OI b, is obtained 
from 



dN olb 



2n 



-dt 



f2 or b dEu 



(2.37) 



^ Stat 

where Etot = E OI b + Bstar is the total energy of the system, 
including both the orbit and the stellar mode. Comparing to 
the fiducial situation in which the mode energy is neglected 
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and dN^l = (f2 or b /27r) (dE OI t, /-Etot ) , we find that the change 
in the number of orbital cycles due to the (a, m)-resonance 
is (Lai 1994) 



(AN, 



orb Ja.m 



t D AE a 



to 



\E a 



(2.38) 



where tzj = |-E or b/-Etot| — D/\D\ is the orbital decay 
time, t or b = 2-k/Q. oy „ is the orbital period, and |-E or b = 
GMM' /(2D) is the orbital energy. The subscript a im- 
plies that to, torb, and |-E or b| are evaluated at the reso- 
nance radius D a . Parameterizing by the resonance index 
(a, m) = (jk, m) and substituting in the expressions for 
AE a , m , we obtain, for the mode with u a > 0, 

5/2 

: i / II \ I i ; <~ \ 

(AWorhJa.m = 



5 / R \ ( Rc 2 
"64^ \GM 2 ) 



(1 + 9) 



2/3 



\Ul a J 



\GM 

7/3 



AE a 



(2.39) 



25 



(Rc 2 Y 1 
6144 \GM J qe a 



E 



(l + g )(2i-l)/6 



( 



(4!-ll)/6 



(2.40) 



If we keep only the dominant Z-term and label the resonance 
by (a, Im), we arrive at 

5 



25 



( Rc 2 



6144 \GM 

(4!-ll)/3 



g(l+ 9 )(2;-l)/3 



X f(^) Wl m \DtS\Q a>lk ?. (2.41) 

& Ql \ III / 

Note that for ui a > 0, the mode energy is always positive, 
and thus (AiV or b) a ,im is negative, which implies that the 
resonance causes the orbit to decay faster. 

Similarly, for the mode with uj a < 0, we find 



(AiY or b)c«,im — + 



25 / Rc 2 
6144 \GM J q(l + g )(2i-i)/3 



\Wg\ 

m 



(4,1- 11)/ 3 



wf m \v 



CO |2| 



(2.42) 



Note that when ui a < 0, the mode energy can be either 
positive or negative (AE a oc uj a e a ), and thus (AN or i,) at i m 
can have both signs. This implies that the resonance can 
either slow down or speed up the inspiral. 

The dominant quadrupolar gravitational wave emitted 
by the binary has frequency f gw = 2z/ rb = £l or b/n. Thus at 
the (a, m)-resonance, the gravitational wave frequency is 

^gw = — K|, (2.43) 
m 

where v a = u a /(2n) is the mode frequency in the inertial 
frame. In the following (§3), we shall express our results in 
terms of f gw instead of \v a \. 

The phase error in the gravitational wave due to the 
resonance is (A$gw)a,m = 47r(A.ZVorb)c«,m- We expect that, 
for (A$g W )a,m ~ 1, the phase error will not affect gravita- 
tional wave detection but could affect detailed analyses of 
the signals, such as parameter extraction (E. Flanagan 1998; 
S. Hughes 1999, private communication). 



The orbital change, or phase error, does not occur in- 
stantaneously at D — D a . The duration of the resonance is 
given by 



(St) a , 



ate 11 ~ late 2 01 b 



2, \ 1/2 ,_2_ 
1 3m 



tDtoib 



1/2 



(2.44) 



This yields the number of orbital cycles during which reso- 
nance occurs, 

f^orb 



(5N OIO )ct,rn — ^ 27T- b ) <k ^°'' 



15 M~ 



5/6 / 1 + g 



1/6 



where M1.4 is the NS mass in units of 1.4M© and recalling 
that f gw = (2/m)u a . The resonance can be considered in- 
stantaneous if (5N OI b) a . m is much less than dN OTO /dlni/ svl , 
the characteristic number of orbital cycles the binary spends 
near orbital frequency f or b = f gw /2, as given by 

1/3 



dN, 



orb 



5c ; ' 



din v e 



1+q 



1 



1927TG 5 / 3 M 5 /3 \ q' A J (TTI/gw) 5 / 3 

229M -,! /3 (^A\ 1/3 ( 



\ 100Hz 



5/3 

1 . (2.46) 



3 ALIGNED OR ANTI- ALIGNED SPIN-ORBIT 

(j3 = 0° OR 180°) 

We now apply the general equations derived in §2 to cal- 
culate the effects of resonance for different types of modes 
in coalescing binary neutron stars. In this section we con- 
sider the cases when fl a and l~2 rb are either aligned (f3 = 
0°) or anti-aligned (j3 = 180°). Various mode properties 
(u) a , Qu,ik, and £ a ) and our procedure to calculate them 
are summarized. The case of general f3 will be discussed in 
§4- 

For (3 = 0° or 180°, we clearly have m = k, i.e., a given 
mode a — {jk} can only be excited once, at the orbital 
frequency given by fcfiorb = 

3.1 F-modes 

For nonrotating NSs, the f-modes have frequencies of order 
(GM/R 3 ) 1/2 , which are too large for any resonance to occur 
during binary inspiral. Rotation can reduce the mode fre- 
quencies, making resonant excitations possible. In this pa- 
per, we shall consider only the f-modes with j = k. Since 
ui a = uj a r ^ — kQ s , these j = k modes are more strongly af- 
fected by rotation. Only the modes with uj a r ^ > need to 
be considered, since rotation actually increases \u> a \ for the 
oj£ r ' < modes, making their resonances impossible (see 
§2.2.1). 

To calculate the f-mode properties relevant for tidal res- 
onance, we shall model the NS as an incompressible Maclau- 
rin spheroid, whose equatorial radius and polar radius are 
ai and 03, respectively. The spin frequency Q s is related to 
the eccentricity e = (1 — a%/a\) 1 ^ 2 of the spheroid via the 
expression 
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l-e 



2\l/2 



(3- 



3e 1 



(3.1) 



where p = 3M/(4tvR 3 ) is the (uniform) density of the 
spheroid, and R = (a^az) 1 ^ is the mean radius. The fre- 
quency of the j = k mode in the rotating frame has been 
obtained by Bryan (1889) and Comins (1979): 

2ke 2 R k 1 



Id, 

77 



M 



= 1±M- 



[(3 - 2e 2 



3e(l 



=2^1/21 



(3.2) 



where the ± sign refers to modes which propagate opposite 
to Q s and in the same direction as fl 3 , respectively, and 



(1 - e 2 ) 1/2 (2k - 1)!! (2p - 2)!! 



2(p-fe) 



(2fc)! 



+ 



(l-e 2 ) 



^(2p-l)!! 

p—k+l 



(1 - e 2 ) 1 ^ 



(3.3) 



For j — k — 2, a more compact expression for the mode 
frequency is given by Chandrasekhar (1969) as 



= n s ± v^TrGpBn - fi 2 



(3.4) 



where Bn is the index symbol (a dimensionless function of 
e) as defined in Chandrasekhar (1969). Also note that for 
Q„ = 0, the mode frequency becomes 



2GM j(j - 1) 
R 3 2j + 1 ' 



(3.5) 



This is essentially the Kelvin mode. As discu ssed above, we 
shall only consider the mode with + sign in (3.2). 
Figure [j] shows the normalized spin frequency, 



v s = v 3 M~l /2 R% 2 = 2170 Cl s Hz, 



(3.6) 



where Rio is the NS radius in units of 10 km and Q s = 
fi s (-R 3 /GM) 1//2 , plotted against the normalized gravita- 
tional wave frequency, 



— A//" -1 / 2 r>3/2 



2170 



(—) \& a \ Hz, (3.7) 



at which resonance occurs [see eq. (2.14)] for several different 
f-modes (j = k = m = 2, 3, 4, 5). Higher order modes (k > 5) 
are not considered since, as shown below, their contribution 
to the orbital change is negligible. The spin frequency at 
which the star becomes secularly unstable (in the absence of 
viscosity) is where ui a becomes negative; for larger fc-values, 
secular instability occurs at lower Q s . Resonance can also 
occur at these rotation frequencies. For completeness, we 
shall include these unstable modes, even though they would 
require very high spin rates. 

The two vertical dashed lines in Fig. ^ indicate the 
frequency range of gravitational wave detectors, i^ gw = 
10 — 1000 Hz. Note that the binary merger occurs at about 
r m in = 3R (for equal mass NSs), corresponding to the max- 
imum gravitational wave frequency 

3/2 



(max) 



1/2 p -3/2 



1181 M^R 



(— y 

V r min / 



Hz. 



(3.8) 



Near this frequency, the distortion to the NS(s) is significant. 
Clearly for resonance to occur before merger would require 



(max) 
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To calculate the tidal coupling coefficient Q a ,ik, we need 
to determine the Lagrangian displacement, £ Q . We begin 
with the Euler equation in the corotating frame for the in- 
viscid fluid: 

9v 

~m 

where 



+ (v • V) v + 2tl s x v = -V*, 



(3.9) 



* = P - + <£-\n 2 s (x 2 +y 2 ). 
P 2 



(3.10) 



$ is the self-gravitational potential. Linearizing and assum- 



ing small perturbations v, $ oc e 1 ' 



,M 



v + 2fi s x v 



-V*. 



eq. (3.9) reduces to 
(3.11) 



Combining this with the continuity equation for an incom- 
pressible fluid, V • v = 0, we find 



d 2 ^ 

dx 2 



a 2 * 
dy 2 



+ i 



d 2 <b 

dz 2 



= 0. 



(3.12) 



A specific set of solutions to this equation, corresponding to 
the j — k mode, is 



A k (x + iy) k V" 



(3.13) 



where Ak is a normalization constant. The Lagrangian dis- 
placement £ = v/(iuig ) can be obtained by substituting 
(j3~il|) into dTfTJ), giving 



kA k 



,(r) (l M 



2Q S 



(x + iy) k 1 (x + iy) e %u 



(3.14) 



Adopting the normalization condition (2.£), we can write 
the displacement vector £ as 

1/2 



(2fc + l)!! 



(2/c)!! A-Kpka\ k a 3 
x (x + ij/) fc_1 (x + iy) e u 



(2k + 1)!! 



1 1/2 



(2k)\\ 4TT P ka 2k a 3 



p p A p J-^ 



Yk k (8AY' 



(3.15) 

(3.16) 
(3.17) 



in Cartesi an, c ylindrical, and spherical coordinates, respec- 
tively. In (3.17), the radial and horizontal wavefunctions are 
given by [neglecting the factor ( — l) fc ] 

1/2 



k 



pa 2k a-j, 



6 



and the spherical harmonic function Ykk is given by 
(-l) fc ,/(2fc + l)< 



Ykk(6,cj>) 



2 k k\ 



47T 



(sin^) e 



k ik$> 



(3.18) 



(3.19) 



Applying appropriate boundary conditions (Comins 1979) 
for the eigenfunction then yields the expression for the mode 
frequency oj^ as given by ([3.2]). 

Using the spherical representation of £ for the f-modes, 
the tidal co uplin g coefficient (for a = {kk} and I = m' = k), 
defined by (|2.12j), can be written as 
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pa\ k a-i 



1/2 



,3 2k-2 

a x pr 



k\Y kk \ 2 + ± 



dO 



dY kk 



(3.20) 



Note that for a given 8, the radial variable ranges from to 
r m ax(#)i witn 

/2\ -1/2 



2 2 \ _1 / 2 

sin cos 



ffll 1 + 



(3.21) 



,V 2 



where // = cos# and £ = (1 - 
then be evaluated explicitly, giving 

Q kk , kk = (g) (J) , 



je. Equation (|3.2CD can 



(3.22) 



recalling the normalization of Q a< i k in units such that G 
M = R = 1. Furthermore, 



T a = -(fc- i)n a 

e Q = a> Q + (fc — l)n s = ii^ r) — fi 3 



(3.23) 
(3.24) 



Substi tutin g equations (3.2), (3.22), and (3.24) into 



(2.41) or (2.42), the change in the number of orbital cy- 
cles due to the resonance can be readily calculated. Figure ^ 
depicts (AN OI h) kk , kk as a function of the normalized gravi- 

For m — k — 2, we have 



(AiVorb) 22,22 



180' 



' 0.895 \ 
> e a ) 



0.69 



( IOOH2 



(3.25) 



where we have normalized i a and Q22,22 to their respective 
nonrotating values (rotation can change these by a factor of 
~ 20 and ~ 1.5, respectively, for the range of Cl s <^ 0.57 con- 
sidered in this paper). The upper sign (negative AiVorb) in 
( 3.25 ) applies when uo a > and j3 = 180°, while the lower 
sign (positive AiVorb) applies when uo a < and /3 = 0°. 
Note that (AiVorb)**,**; diverges for k — 2 as v SV! decreases. 
The reason for this is that at larger orbital radii (small reso- 
nant frequencies), the orbit decays more slowly, and the bi- 
nary spends a longer time in resonance, giving rise to larger 
AiV or b • For higher fc-resonances, the same effect is also oper- 
ative, but the tidal potential is weaker so that the net effect 
is that (AiVorb)**, ** oc i/gt* u '/ 3 decreases with decreasing 



into account when constructing theoretical wave templates 
in the search for gravitational wave s fro m the se sy stems. 

It can also be seen from eqs. ( |3.l[ ) and (3.2) that the 
NS spin frequency and mode frequencies depend on the ra- 
dius of the NS and thus the NS equation of state (EOS). 
Current equations of state allow NSs with R = 15 km (e.g., 
Shapiro & Teukolsky 1983). This implies that the spin and 
mode frequencies are lowered, and resonance occurs at these 
lower frequencies as illustrated by Figure ^. Figure ^ shows 
the corresponding (AN OI h) kk , kk as a function of the v svr 
for the R — 15 km NS. The lower limit to the spin fre- 
quency required for resonance is set by the associated grav- 
itational wave frequency given by eq. (3.S). This limit is 
indicated by the dotted lines in Fig. |^and Fig. ^ and yields 
v s <; 330 Hz for the (22, 2)-resonance, v a ^ 260 Hz for the 
(33, 3)-resonance, u s <; 200 Hz for the (44, 4)-resonance, and 
v s <; 150 Hz for the (55, 5)-resonance. 

It is necessary to check that the higher-i contribu- 
tion to the (kk, m = fc)-resonance is negligible. Except 
for the Q s = case, the I = k + 2, k + 4, • • ■ tides can 
also contribute to (AiV rb)*fc,fc- The mode frequency is the 
same, but the new tidal coupling coefficient must be cal- 
culated. We shall consider just the (22,42)-resonance and 
the (33,53)-resonance here since the higher or der t erms are 
much smaller. Using the mode eigenfunction (3.17), we find 



Q'2 



45 
112 
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2tt (1 - e 2 



1/2 



; 1 -5 M ' 4 + 6 M ' 2 - 1 , , 



(3.26) 



and 



35 
96 



11 



-, 1/2 



7T (1 



7/i' 6 - 15m' 4 + V 2 



■ dfj! ' . 



(3.27) 



Both expressions can be evaluated numerically. In Fig. |^ the 
results for the orbital change (AiV or b)22,42 is plotted, indi- 
cating that it is indeed much smaller than (AiV or b)22,22[J} 
The orbital change (AiV or b)33,53 is not shown since it is be- 
low the limit of interest for this paper. 

3.2 G-modes 



It is evident from Fig. g that the (jk, m) = (22, 2) 
and (33, 3)-resonances contribute significantly to the or- 
bital phase evolution. But, to attain such resonances dur- 
ing the binary inspiral, the NS of canonical mass (IAMq) 
and radius (10 km) must be rotating very rapidly, i.e., 
v s 'Z 710 Hz for the (22, 2)-resonance, v s ^ 570 Hz for the 
(33, 3)-resonance, v s Jj 460 Hz for the (44, 4)-resonance, and 
v s <; 380 Hz for the (55, 5)-resonance (see Fig. |l]). Such 
rapid rotation is physically possible but may be astrophys- 
ically unlikely. Observations of NS binary systems and the 
current theories regarding the formation of these systems 
indicate that such rapid rotation of the NS is not possible 
(see §5). If such resonances indeed occur, then their effect 
on the orbital phase of the binary inspiral must be taken 



Gravity modes in neutron stars arise from the composi- 
tion (proton to neutron ratio) gradient in the stellar core 
(Reisenegger & Goldreich 1992; Lai 1994), density disconti- 
nuities in the crust (Finn 1987) , as well as thermal buoyancy 
associated with finite temperatures, either due to internal 
heat (McDermott et al. 1988) or due to accretion (Bildsten 
& Cutler 1995). In this paper, we shall focus on the core 



II Note that to calculate the total orbital change due to the 
(jk,m) = (22, 2)-resonance, one needs to use eq. (2.40) which 



includes contributions from all possible i-tcrms. But the smallncss 
of (A7V orb )22,42 implies (A7V orb )22,2 =i (A7V orb )22,22 to a good 
approximation. The quantity (A7V or b)22,42 will also be used in §4 
as a fiducial number for the (22, 44)-resonance. 
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g-modes since the gradient of proton to neutron ratio in the 
NS interior provides the largest buoyancy. For a nonrotating 
NS, the typical frequencies of low order g-modes are around 
100 Hz and are smaller for higher order modes. Thus we ex- 
pect that even modest stellar rotation rate (~ 100 Hz) can 
change the mode properties dramatically. 

The core g-modes of rotating NS have been calculated 
by Lai (1998) using the traditional approximation. Two NS 
models are considered based on different EOS. Model UU 
gives a nonrotating NS radius R — 13.47 km for a M = 
IAMq, and its lowest order (n = 1) j = 2 g-mode has 
frequency v a = 148 Hz and Q 2 2,22 ^ 0.62 x 10~ 3 . Model 
AU gives a nonrotating NS radius R — 12.37 km for a M = 
1.4Mq, and its lowest order (n = 1) j = 2 g-mode has 
frequency v a = 72 Hz and Q22,22 — 0.11 x 10 -3 . 

Figure ^ shows the spin frequency, v s , as a func- 
tion of the observed gravitational wave frequency, f gw = 
(2/m)\u a \ = \v a \, for the m = 2 resonance of the {jk} = 
{22} g-modes (n = 1). The uj a r ' > > mode propagates in 
opp osite direction as the spin [see the discussion following 
eq. ( 2.20 )]. As v 3 increases, the mode frequency, u a , in the 
inertial frame decreases; v a can even become negative when 
u s exceeds a certain critical value (95 Hz for model UU and 
47 Hz for model AU). This corresponds to the onset of sec- 
ular instability for the g-mode (Lai 1998) . We include these 
unstable modes because, as shown by Lai (1998), viscosity 
tends to stabilize these modes for typical parameters of the 
NS. 

Figure |^ also shows the results for the spin-prograde 
mod e, w hich has LUcP < [see the discussion following 
eq. (2.20)]. These are calculated in a similar manner as in 



Lai (1998). These modes always have negative ui a , but they 
have positive energy and thus are secularly stable. Rotation 
clearly increases the gravitational wave frequency at which 
resonance occurs for these modes. 

Figure [] shows the orbital change, (AA r or b)22,22, due to 
the g-mode resonance as a function of the observed gravi- 
tational wave frequency, f gw . The mode parameters Q22.22 
and e a ha ve be en o btain ed using the method of Lai (1998). 
Using eq. (2.41) or (2.42) (depending on the sign of we 
find 



(AJV„. 



180° 




=F0.0086 



?(! + «) 



/22.22 



6.2 x 10" 



50Hz 



(3.28) 



for model UU and 
/l80 c 



(AiVo 



V o° 



=F0.00045 



9(1 + 3) 



1.1 x 10- 4 



50Hz 



(3.29) 



for model AU. We have normalized e a and Q22,22 to their 
respective nonrotating values (rotation can change these by 
a factor of ~ 2 for the range of C l s <^ 0. 3 con sidered in this 
paper). The upper sign in (3.28) and (3.29) applies when 
to a > and (3 = 180° 



while the lower sign applies when 
io a < and (3 = 0°. However, when w£ < and u a < 0, 
e a < 0; therefore AJV or b < (a decrease in the number of 
orbits) . 



The large dots in Fig. ^| correspond to the case when 
the NS is nonrotating. We see that retrograde rotation (f3 = 
180°) can increase |(AjV or b)22,22| for the spin-retrograde 
mode by shifting the resonance to smaller orbital frequency, 
where the orbital decay is slower. For the spin-prograde 
modes, rotation decreases |(AJV orb )22,22| from the nonrotat- 
ing value. It is clear from Fig. ^| that the orbital change is 
small, especially in comparison to the orbital change gener- 
ated by f-modes. However, the g-mode resonances can occur 
more easily than the f-modes since the g-mode resonance 
does not require the rapid rotation of the NS. 

The above r esults apply for particular neutron star 
models. As in eq. (3.25), we have the scaling: (AA r or b)22,22 oc 



M 1 4 Rfo; thus the orbital change is larger for low-mass 
NSs which typically have larger R/M. 



3.3 R-modes 

For nonrotating stars, r-modes are "degenerate" toroidal 
modes with zero frequency and zero density perturbation. 
Rotation breaks the degeneracy via the Coriolis force, giv- 
ing rise to finite-frequency r-modes. For the analysis of the 
r-modes in this section and in §4.3, we shall consider the 
slowly-rotating [Cl s = n s /(GM/R 3 ) 1/2 < 1] incompressible 
NS models. As shown by Provost et al. (1981), only the j = k 
r-modes exist for barotropic stars. 

Following Saio (1982) and Kokkotas & Stergioulas 
(1998), it is convenient to introduce a new radial coordi- 
nate a which is related to the radial distance r through 
r — a(l + e), where eia.,6) oc Cl 2 s represents the centrifu- 
gal distortion of the star. In the new coordinates, all unper- 
turbed quantities are functions of a only. The displacement 
eigenfunction of the a — {kk} r-mode can be written as 



¥ - [ Q > Kkk iLik~ Kkk ire 



— I Sk+ik, H k+lk — , H k+ i k -^—j Yfe+ifc, (3.30) 

where the first term represents the toroidal (axial) compo- 
nent and the second term represents the spheroidal (po- 
lar) component. The r-mode is predominantly toroidal, and 
the radial eigenfunctions scale with the rotation rate as 
Kkk ~ Cls, while Sfc+ifc, H k +ik ~ ^a- For a uniform den- 
sity star, the mode frequency in the rotating frame is given, 
to order f2 3 , by (Saio 1982; Kokkotas & Stergioulas 1998) 



2 fi. + 5 * ( * +1 £- 8 n; 



fc + l"" 1 (fc + 1) 4 

In the inertial frame, this gives 



(fc + 2)(fc- 



k- 



(k- 



(3.31) 



(3.32) 



Note that the r-mode always has positive u„ and nega- 
tive Lo a , i.e, the wave pattern propagates opposite to spin 
in the rotating frame while it becomes spin-prograde in the 
inertial frame, indicating that the r-mode is secularly un- 
stable, although viscosity tends to stabilize the mode for 
small Cl s (Andersson 1998; Andersson et al. 1998; Friedman 
& Morsink 1998; Lindblom et al. 1998). 

The eigenfunction (3.30) implies the density perturba- 
tion for the r-mode Sptk oc Yk+ik- Thus the tidal coupling 
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coefficient, Qkk,lk, is nonzero only for I = k + 1. To evaluate 
Qkk.lk, we write the pressure perturbation of the r-mode as 



5pkk = pga( k +ik(a)Y k +ik, 



(3.33) 



where g is the gravitational acceleration. The dimensionless 
function £k+ik(a) is given by (Kokkotas & Stergioulas 1998) 



Cfc+ifc(a) = (f 3 (~j 

and the corresponding radial eigenfunction is 
.(k + l)V2fc + 3 fa^ k - 1 



Kkk(a) = i- 



4A- 



(1) 



The Eulerian density perturbation is given by 
Spkk = — ^fefc ■ Vp = apSk+ikY k +ik5(a — R). 



(3.34) 



(3.35) 



(3.36) 



But the boundary condition at the stellar surface gives 
Sk+ik(R) = Cfc+ifc(-R)> which is obtained by requiring the 
Lagrangian pressure perturbation to vanish at the surface. 
We then find 



Q 



kk,k + 



ik = pfl s R 



fc+4 



(3.37) 



Note that t his v alue corresponds to the normalization as 
defi ned by ( 3.34 ) . To obta in the proper normalization [see 
eq. (2.9)], we use eq. (3.3C) and find, to leading order in Cl s , 



d 3 x P i- e 



k(k + l)p / daa 4 \Kkk{a)\ 2 
Jo 

(k + 1) 5 



16k 



- P R° 



(3.38) 



Thus the renormalized Qkk,k+ik, expressed in units where 
G = M = R= 1, is 



Q 



fcfc,fc + lfc 



12k 



n(k + l) b 



1/2 



ni 



(3.39) 



Thus tidal co uplin g will be weak for slow rotation. With the 
£a given by ( 3.3C ), we ca n sim ilarly calculate the function 
T a and e a [see eqs. (2.11) and (2.3C)], yielding 



£ a — LU a + kti s = (2>£\ 



(3.40) 
(3.41) 



The next order correction to T a is of order Cl 5 , and can be 
neglected. 

Now consider resonant excitation of the a = {kk} r- 
mode. For aligned spin-orbit we require m = k, while to have 
non- vanishing Qkk,ik we require I — k+1. Since Wk+ik = 0, 
we conclude that for (5 — or 180° , there is no effective tidal 
excitation of the r-mode (at least for the barotropic stellar 
model considered here) . Resonant excitations of r- modes can 
only occur for misaligned spin-orbit configurations, which we 
consider in §4.3, using the properties of the r-mode derived 
above. 



4.1 F-modes 

We again focus on the j = k modes. For a given a = {kk} 
mode, in addition to the (a, m) = (kk, fc)-resonance consid- 
ered in §3.1, many new resonances become possible for gen- 
eral (3. This c an be seen by inspecting the resonance con- 
ditions ( 2.14)-(2.16): For Qkk,ik to be nonzero, we require 
/ = k, k + 2, k + 4, • • • To have nonzero tidal potential 
{Wim 7^ 0) we must then have m — k, fc±2, k+A, ■ ■ ■ (m > 0). 
For example, the j = k = 3 mode can be resonantly ex- 
cited at orbital frequencies Q OI b = |^q|, ^orb = |w a |/3, and 
f2 orb = \cu a \/5, etc. 

The properties of the f-modes have already been calcu- 
lated in §3.1. For uj a > (the spin-retrograde mode), the 
relevant Wigner ©-functions are 



V 



(fc) 



K) 



v 



fc,-fc+2 



'y/k(2k- 



v 



fc,-fc+4 



lka /*(2fc-l)(Jfc-l)(2fc-3) 



X ( C ° S f) ( Sin f 
Vl% 2) =e ika (sinf) 2fc [(sinf) 4 -2(2fc + 2) 

x (cos0 2 (sin^) 2 + (* + 1)(2* + 1) (cos!)' 



(4.1) 
(4.2) 

(4.3) 



(4.4) 



From eq. (2.41), we find that the orbital change due to the 
(kk, Zm)-rcsonance at inclination angle /3 can be written as 

(AN orb ) kk , lm (f3) = (AiY orb )fcfc, !fc (180°) x f k k,im(J3), (4.5) 

where (AiV or b)fcfc,zfc(180 o ) is the value for anti-aligned spin- 
orbit, and the fkk,lm is a dimensionless function 



,(4i-ll)/3 



Wi-kJ <~ k ~ ml KmJ (4 ' 6) 

Figure ^ shows f k k,lm(fi) as a function of the inclina- 
tion angle for a number of different resonances. This fig- 
ure should be read together with Fig. |^ to determine the 
actual orbital change at different v gv and f3. The results 
for the (33,53) and (33,55)-resonances are not shown since 
their contributions to the orbital change are below the lim- 
its of interest. It is evident from Fig. M that even when the 
rotation of the NS is mostly prograde with respect to the 
orbital motion ((3 < 90°), there still exist modes that can 
be resonantly excited and thus contribute to orbital energy 
dissipation. Also of note is the large orbital change caused 
by the (55,51)-resonance relative to the (55,55) and (55,53)- 
resonances. It is thus possible for this resonance to produce 
an orbital change which is comparable to those produced by 
the (44,4m)-resonances. 

For spin-prograde modes (u) a < 0), the above results 
still apply except that one should replace (AA r or b)fcfc,(fc(180°) 



4 NON-ALIGNED SPIN-ORBIT 

In this section we consider tidal resonances for general spin- 
orbit inclination angle, j3. 



** This is because the eigenfunction of the a = {kk} mode is 
even (odd) with respect to the stellar equator for even (odd) A;; 
the tidal potential must have the same symmetry in order to 
produce nonzero coupling. 



© 0000 RAS, MNRAS 000, §-0 



Resonant Tidal Excitations of Rotating Neutron Stars 



11 



by (AiV or b)fefc,(fc(0 o ) and, in the expression for fkk,lm(P), re- 
place P by (180° -/?). 



4.2 G-modes 

The g-mode resonance considered in §3.2 corresponds to the 
(a,lm) = (22,22). The function fkk,im(P) derived in §3.1 
equally applies to the g-mode resonance and is plotted in 
Fig. M. Although other resonances of g-modes are also in- 
duced, their effects are below the limits of interest. 



4.3 R-modes 

Recall that only the j — k r-modes exist for barotropic stars, 
and Q kk ,ik 7^0 only if I = k + 1, in which case Q k k,k+ik 
is given by (3.39). Resonances at m are restricted to m = 
k + 1, k — 1, k — 3, . . . (m > 0) since Wk+im = otherwise. 
We consider only the (22,3m)-resonances (for m = 3, 1) since 
the orbital effect decreases with increasing I. The relevant 
Wigncr ©-functions are 



7,(3) 
£>(3) 



V6s 



P 



( P\ 
cos — 

V 2) 



10 



( P\ s L/ . P\ 2 ( py 

cos— 2 sin- — cos — 
V 2) \ 2) \ 2) 



(4.7) 



(4.8) 



For the (22,33)-resonance, we find, to leading order in 
Cls, that the gravitational wave frequency is given by 



8 

— Vs 

9 



(4.9) 



where in th e second equality, we have dropped the Cl^ term 
in eq.(3.32). The orbital change is given by 



(AiVo: 



3.0 x 10 



_ 4 ( Rc 2 ' 

\GM j g(l + g)s/3 

P\ 2 ( P" w 



- ('■) 



*(-f)(-IJ 



+ 4.6 x 10 -5 Mi 



20/3 D 10 



10" 



100Hz 



g(l + g) 5 / 3 

10/3/ fl\2/ o 10 

(sinf ) (cosf ) . (4.10) 

Recall that the r-mode has negative energy, and as a result, 
(AN OI b) a ,im is positive, i.e., the resonance slows down the 
inspiral. At the maximum j/ gw = 1000 Hz (or v s ~ 1125 Hz) 
and P = 48.2° (which would give the maximum orbital 
change), we have 



(AJVort,)^ ~ 0.010 M-f /3 Rl° 



9(1 + ' 



\5/3 ■ 



(4.11) 



The orbital change increases with increasing spin frequency. 
But note that our analysis of r -mod e is based on the as- 
sumption of slow rotation, thus (4.11) should be considered 
as an estimate. 

For the (22,31)-resonance, we have 



= 2\u a 



8 

— fa 

3 



(4.12) 



and 



(AN o: 



+ 4.3 x irr 
P\ 2 ( P^ 6 



sm— cos— 2 sin— — cos— 
V 2) V 2) \ 2) \ 2) 



+ 6.6 x 10" 

P^ 



20/3 D 10 
-Kio 



q(l + q) 5 ^ 
P\ 2 ( P^ 2 



(sin.) (cos-) 2 (sin-) - (cos-) 



100Hz 

2 



10/3 



(4.13) 



Again, at the optimal parameters i/ gw = 1000 Hz (or v a 
375 Hz) and P = 34.4°, we have 



(AAorb)^^ 



2.0 x 10~ 4 M- 4 20/3 < 



)5/3 ■ 



(4.14) 



Thus we conclude that for canonical NS parameters, 
the (22,33) and (22,31)-resonances (and indeed all r-mode 
resonances) produce orbital changes that are below the lim- 
its of interest. However, if we consider a larger NS radius, 
such as R — 15 km for a M = 1.4M Q NS, which is allowed 
for some nuclear equations of state (e.g, Shapiro & Teukol- 
sky 1983), then the orbital change for the (22, 33)-resonance 
may become significant (AiV or b ~ 0.6). 



5 SUMMARY AND DISCUSSION 

In this paper, we have presented a systematic study on the 
resonant mode excitations of rotating neutron stars in coa- 
lescing binaries and their effects on the gravitational wave- 
form. The essential effect of rotation is that it can modify 
the mode frequency in the inertial frame, thereby changing 
the strength of the g-mode resonance, which already exists 
for nonrotating neutron stars, and also making a variety of 
new resonances involving f-modes and r-modes possible. 

We find that for the f-mode resonance to occur during 
the last few minutes of the binary inspiral, with the emitted 
gravitational wave frequency 10 Hz < z/ gw < 1000 Hz (the 
sensitivity band of the interferometric gravitational wave de- 
tector), the neutron star must have very rapid rotation (see 
Fig. hi), e.g., assuming canonical neutron star mass and ra- 
dius, M = 1.4M and R = 10 km, respectively, v s ^ 710 Hz 
for the (jk, m) = (22, 2)-resonance (see the beginning of §2.2 
for notation) and v a ii 570 Hz for the (33, 3)-resonance. Ob- 
servations (e.g., PSR 1913+16 has spin period 59 ms, and 
PSR 1534+12 has 38 ms) and our current understanding of 
the formation of neutron star binaries seem to indicate that 
neutron stars in compact binaries (with another neutron star 
or a black hole as a companion) can never achieve such a 
rapid rotation, although such a rotation rate is physically 
allowed. For a M = 1.4Mq and R = 15 km neutron star, 
the critical frequencies are lowered, e.g., v a ^ 330 Hz for the 
(22, 2)-resonance and v s J; 260 Hz for the (33, 3)-resonance 
(see Fig. |§]). If the neutron star in a coalescing binary in- 
deed has the required rapid rotation, then the induced or- 
bital change due to the f-mode resonance is significant (see 
Fig. ^ and Fig. ^) and must be included in the templates of 
waveforms used for searching gravitational wave signals. 

The resonance of g-modes is strongly affected by even a 
modest rotation (see Fig. ^) . Although there is large uncer- 
tainties associated with the property of the g-mode (since 
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it depends on the symmetry energy of nuclear matter), we 
show that for a canonical IAMq neutron star, the orbital 
change A7V or b due to the g-mode resonance is below the 
10 -2 level (see Fig. |^). Such a small phase error in the wave- 
form is unimportant for the search of gravitational wave sig- 
nals from the binary (E. Flanagan 1998, private communi- 
cation). However, we note that if we consider very low-mass 
(M ^ 0.5Mq) neutron stars, which have large R/M ratio, 
the orbital change due to g-mode resonance may well be 
significant. Of course, although such low-mass neutron stars 
are physically allowed, there is no astrophysical evidence for 
their existence. 

Resonant excitations of r-modes occur for misaligned 
spin-orbit configurations. Since the r-mode always has neg- 
ative energy, the resonance actually transfers energy to the 
orbit, slowing down the inspiral. Since the tidal coupling of 
the r-mode depends strongly on the stellar rotation rate, we 
find that for canonical neutron star mass and radius (IAMq 
and 10 km, respectively), the orbital change AN m b is neg- 
ligible. However, if we consider a larger stellar radius, such 
as 15 km for a IAMq neutron star (as is allowed by some 
nuclear equations of state), then AN OI b may become impor- 
tant, esp ecially in th e hig h frequency band (j/ gw ~ 1000 Hz) 
[see eqs. ( |4.1o| ) and (jUl)], 

Finally, we note that our calculations of the properties 
of the f-modes and r-modes are based on incompressible neu- 
tron star models. Our results for the r-modes should break 
down for very high rotation rates. Since neutron stars are 
not exactly barotropic, other r-modes, in addition to the 
j = k mode considered in this paper, can also exist. It is 
desirable to study the resonance effects using more realistic 
neutron star models. 
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Figure 1. Normalized spin frequency, P s = u s M 1 l^ 2 R^g 2 , of the NS as a function of the normalized gravitational wave frequency, 
P gw , for the f-mode (jk, m)-resonance (taking j = k = 2—5), where P gw is related to the f-mode frequency by P gw = Vgw M 1 ^ 2 R^q 2 = 
(2/m)\v a | M-^ ^ 2 Riq 2 ■ The labels in the figure give the values of (jk, m). The solid lines are for the oia' > and ui a > (spin-retrograde) 
stable modes. The dot-dashed lines are for the oj£ r ' > and u) a < (spin-prograde) unstable modes. Note that the (22,4)-resonance is 
possible only for misaligned spin-orbit configurations (/3 ^ 0°, 180°). It is also clear the unstable modes occur at extremely high spin 
rates. The vertical dashed lines indicate the frequency range of gravitational wave detectors, u gvi = 10-1000 Hz. 



© 0000 RAS, MNRAS 000, 0-|l| 



14 W. C. G. Ho and D. Lai 



Figure 2. Orbital cycle change |(AA r or b)jfc i ; m |, with M = IAMq, R = 10 km, q = 1, as a function of the normalized gravitational wave 
frequency, P gw = i/ gw M 1 \^ 2 R\(^ = (2/m)|^ a | M 1 ^ 2 R^q 2 , for the f-mode (jk, Zm)-resonance. The labels in the figure give the values of 

(r) 

(jk,lin). The solid lines are for the uj„ > and u) a > (spin-retrograde) stable modes with (AA r or t,)jfc > ; m (180°) < 0. The dot-dashed 
lines are for the oj^ > and uj a < (spin-prograde) unstable modes with (A7V or b) :) 7 :i ( m (0 o ) > 0. Note that both the (22,22) and (22,42) 
curves belong to the same (22,2)-resonance. The (22,42) curve serves as a baseline for the (22,44)-resonance (see §4 and Fig. 7). Also 
recall from Fig. 1 that the unstable modes require extremely high spin rates. The large dots with labels indicate the spin rate at that 
gravitational wave frequency. The vertical dashed lines indicate the frequency range of gravitational wave detectors, f gw = 10—1000 Hz. 
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Figure 3. Spin frequency, u s , of the R = 15 km NS as a function of the gravitational wave frequency, v gw , for the f-mode (jk,m)- 
resonance (taking j = k = 2-5), where u gvJ is related to the f-mode frequency by i> gw = (2/m)\u a \. The labels in the figure give the values 
of (jk,m). The solid lines are for the £j„ > and u) a > (spin-retrograde) stable modes. The dot-dashed lines are for the u)^ ] > 
and ui a < (spin-prograde) unstable modes. Note that the (22,4)-resonance is possible only for misaligned spin-orbit configurations 
(/3 ^ 0°, 180°). It is also clear the unstable modes occur at extremely high spin rates. The vertical dotted line at i^ gw = 640 Hz denotes 
the maximum gravitational wave frequency allowed by eq. (3.8), which sets a lower limit to the spin frequency. The vertical dashed lines 
indicate the frequency range of gravitational wave detectors, u gvl = 10—1000 Hz. 
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Figure 4. Orbital cycle change \(AN OI b)jk^ rn \, with M = 1.4Mq, R = 15 km, q = 1, as a function of the gravitational wave frequency, 
fgw = (2/m)\v a \, for the f-mode (jk, im)-resonance. The labels in the figure give the values of (jk,lm). The solid lines are for the 
> and ui a > (spin-retrograde) stable modes with (AAf or b)jfc > i m (180°) < 0. The dot-dashed lines are for the oj^ r ' > and 
co a < (spin-prograde) unstable modes with (AAf or ] :i ) :f 7 S! i m (0 o ) > 0. Note that both the (22,22) and (22,42) curves belong to the same 
(22,2)-resonance. The (22,42) curve serves as a baseline for the (22,44)-resonance (see §4 and Fig. 7). Also recall from Fig. 2 that the 
unstable modes require extremely high spin rates. The vertical dotted line at v gw = 640 Hz denotes the maximum gravitational wave 
frequency allowed by eq. (3.8). The large dots with labels indicate the spin rate at that gravitational wave frequency. The vertical dashed 
lines indicate the frequency range of gravitational wave detectors, f gw = 10-1000 Hz. 
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Figure 5. Spin frequency, u B , of the NS as a function of the observed gravitational wave frequency, v gw , for the (jk,m) = (22,2)- 
resonance, where u gv , is related to the g-mode frequency, u a , by u gv = (2/m)\v a \ = \v<x\- UU and AU are two models of the NS EOS 
(see text for the NS parameters) given in Lai (1998). The solid lines are for the ' > and ui a > (spin-retrograde) stable modes. 

(r) z \ (r) 

The dot-dashed lines are for the u) a > and (jj a < (spin-prograde) unstable modes. The dotted lines arc for the uj^ < and ui a < 
(spin-prograde) stable modes. The vertical dashed line indicates the lower limit of the frequency range of gravitational wave detectors, 
u gv = 10-1000 Hz. 
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Figure 6. Orbital cycle change |(A.iV or b)22,22 |i with q = 1, as a function of the observed gravitational wave frequency, f gw = (2/m)\u a \ = 
\u a \, for the g-mode (jk,lm) = (22, 22)-resonance. UU and AU are two models of the NS EOS (see text for the NS parameters) given 
in Lai (1998). The solid lines are for the lu^, > and u) a > (spin-retrograde) stable modes with (A7V or t,)22,22 (180°) < 0. The 

( r) 

dot-dashed lines are for the ui^, > and ui a < (spin-prograde) unstable modes with (AAf or i D )22,22 (0°) > (Note that the result 

( r) 

for model AU is below the limit of the figure). The dotted lines are for the ui^, < and u) a < (spin-prograde) stable modes with 
(AA r or b)22,22(0°) < 0. The large dots indicate where the NS is nonrotating. The vertical dashed line indicates the lower limit of the 
frequency range of gravitational wave detectors, u gvl = 10—1000 Hz. 



© 0000 RAS, MNRAS 000, §-0 



Resonant Tidal Excitations of Rotating Neutron Stars ... 19 



Figure 7. fkklmifi) as a function of spin-orbit inclination angle, /3. f^k im{P) is related to the induced orbital cycle change, 
(AiV orb ) fefc j m (/3), by (AAf orb ) fcfc , im (/3) = (AAf orb ) fefc j k (180°) X f k k,hnW)- The labels in the figure give the values of (jk,lm). 
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